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We present branching fraction measurements for the radiative decays B"*" —^ p'^'y, B° — > p°7, 
and B" —* ijjy. The analysis is based on a data sample of 465 million BB events collected with 

the B/lRlR, detector at the PEP-II asymmetric-energy B Factory located at the Stanford Linear 
Accelerator Center (SLAC). We find B{B+ p+7) = (1.20t^;37 ± 0-20) x 10"^ B{B° p^-j) = 
(0.97lo:22 ± 0.06) X 10"®, and a 90% C.L. upper limit 0^7) < 0.9 x 10"^ where the first 

error is statistical and the second is systematic. We also measure the isospin-violating quantity 

r(s+ ^ p+7)/2r(s° ^ p°7) - 1 = -oA3l:°oil ± o.io. 
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I. INTRODUCTION 

Within the standard model (SM), the radiative de- 
cays P~^l, P^l, and Luj ^ proceed 
mainly through a, b ^ d'y electroweak penguin ampli- 
tude with a virtual top quark in the loop. Hence, the 
decay rates depend on the magnitude of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Vtd- The 
branching fraction results from recent next-to-leading or- 
der calculations are listed in Table I. While these ex- 
clusive decay rates have a large theoretical uncertainty 
dominated by the imprecise knowledge of the form fac- 
tors, some of this uncertainty cancels in the ratio of 
B p{(jj)^ to B K*'y branching fractions. This ra- 
tio provides a constraint on the ratio of the CKM matrix 
elements |V(d/Vts|, which can also be obtained from the 
ratio of Bd and Bg mixing frequencies [4]. Physics be- 
yond the SM could affect differently B p{^)l and 
Bd/ Bs mixing, and hence create inconsistencies between 
the results obtained from the two methods. 

The ratio of B ^ p{^)l to i? — !■ K*^ branching frac- 
tions is related to |Vf(j/Vis| [5] via 



TABLE I: Recent predictions of the branching fractions. 



B[B ^ p(a;)7] 
B{B K*-f) 



S 



Vtd 



Vts 




(1) 

The coefficient 5 is 1 for /9+ and | for p° or w, m is the 
particle mass, Cp(w) is the ratio of the form factors for the 
decays B p{co)^ and B K*j, and ARp(^^^ accounts 
for differences in decay dynamics, including weak annihi- 
lation contributions. The precision of the |14d/Vts| deter- 
mination can be improved by using an average branch- 
ing fraction for B /o(w)7 decays. Within the SM, the 
isospin asymmetry between /9+7 and B° p°7 

is dominated by weak annihilation contributions, and is 
expected to be small; on the other hand, the asymmetry 
between B° — » p'^^ and — * loj can be sizable, due to 
the difference in the form factors [1, 3]. 

We report an updated study of the decays i?+ 
p+7, B° p°7, and S° (jj"f based on 465 million 
BB events, corresponding to an integrated luminosity 
of 423 fb~^, a data sample 25% larger than that use 



Mode 




Branching fraction (xlO ®) 
Ref. [1] Ref. [2] Ref. [3] 


B+ ^ 


p+7 


1.41 ±0.27 1.58iSi| 1.16 ±0.26 




p\ 


0.69 ±0.12 0.76lo;23 0.55 ±0.13 






0.55 ±0.09 0.44 ±0.10 



in our previous publication [6]. In addition, we reduce 
backgrounds considerably by using a multivariate algo- 
rithm based on bootstrap-aggregated (bagged) decision 
trees (BDTs) [7] and additional discriminating variables 
to separate signal from background. 



II. THE BABAR DETECTOR AND DATA SET 

The data sample is collected with the BaBar detec- 
tor at the PEP-II asymmetric-energy e'^e~ storage ring 
at a center of mass (CM) energy near -^i = 10.58 GeV, 
corresponding to the T(4S') resonance (on-resonance). 
Charged particle trajectories and energy loss (dE/d.T) are 
measured with a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH) in a 1.5 T magnetic 
field. Photons and electrons are detected in a CsI(Tl) 
crystal electromagnetic calorimeter (EMC) with photon 
energy resolution as/E = 0.023(i;/GeV)-i/'* ® 0.019. A 
ring-imaging Cherenkov detector based on the detection 
of internally reflected Cherenkov light (DIRC) provides 
information for charged particle identification. The K-n 
separation in the DIRC is above 4c7 at laboratory mo- 
menta up to 3 GeV/c. In order to identify muons, the 
magnetic fiux rctTirn is instrumented with resistive plate 
chambers and limited streamer tubes. A detailed descrip- 
tion of the detector can be foimd elsewhere [8] . 

We use a GEANT4-based [9] Monte Carlo (MC) simu- 
lation to model the BaBar detector response, taking into 
account the varying accelerator and detector conditions. 
Dedicated signal and background MC samples are used 
to optimize selection criteria, to obtain signal efficien- 
cies and to validate the analysis. Data control samples, 
including 41 fb~^ of data collected about 40 MeV below 
the BB production threshold (off-resonance), are used to 
study backgrounds coming from continuum e+e" — > qq, 
with q = u,d, s, c. 
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III. EVENT RECONSTRUCTION AND 
BACKGROUND SUPPRESSION 

The decays B p(w)7 are reconstructed by combin- 
ing a high-energy photon with a vector meson recon- 
structed in the decay modes /9+ 7r+7r°, p° tt+tt", 
and w — > 7r+7r~7r°. The dominant source of background 
is coming from continuum events that contain a high- 
energy photon from 7r° or r] decays or from initial-state 
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radiation (ISR). There are also significant backgrounds 
from B meson decays. Ttie decays B — > K*'y, K* Kn, 
can mimic the signal when the kaon is misidentified as 
a pion. Decays of B — » [p / u:){-k^ / rj) with a high-energy 
photon from the 7r° or rj decay also mimic the signal. 
In addition, there arc other B backgrounds originating 
mainly from B and B X[-k^ /r}) decays. 

The event selection and background suppression are 
performed in two steps. We apply a set of loose selection 
criteria to select well-measured photons and charged pi- 
ons and to reject background events that are kinemat- 
ically very different from the signal events. For events 
that pass the loose event selection criteria, we then use 
the BDT technique to further reduce background. 



A. LOOSE SELECTION 

We reduce background contributions from continuum 
processes by considering only events for which the ratio 
i?2 of sccond-to-zcroth Fox- Wolfram moments [10], cal- 
culated using the momenta of all charged and neutral 
particles in the event, is less than 0.7. 

A photon candidate is identified cluster of energy 
deposited in contiguous EMC crystals, and not associ- 
ated with any charged track. The high energy photon 
must have energy 1.5 < < 4.4 GeV in the labora- 
tory frame and 1.5 < E* < 3.5 GeV in the CM frame, 
be well-contained within the EMC acceptance with polar 
angle —0.74 < cos9 < 0.93, and be isolated by at least 
25 cm at the entrance of the EMC from any other pho- 
ton candidate or charged track. The distribution of the 
energy deposition is required to be consistent with that 
of a photon shower. 

Charged-pion candidates arc selected from well- 
reconstructed tracks that have at least 12 DCH hits used 
in the track fit and a minimum momentum transverse 
to the beam direction of 100 MeV/c. The tracks are re- 
quired to originate near the interaction point (IP): the 
distance of closest approach to the IP must be less than 
10 cm along the beam direction and less than 2 cm in 
the plane perpendicular to the beam direction. The tt^ 
identification is based on a likelihood Li computed for 
particle hypothesis i{= Tr,K,p) using dE/dx measured 
in the SVT and DCH and the information of Cherenkov 
photons detected by the DIRC. The selection criteria are 
optimized to reject charged kaons produced in _B ^ K*j 
decays. The pion candidates in —> toj must have 
Lk/{Lk + L^) < 0.5 and Lp/{Lp + L^) < 0.98 and 
must not be consistent with being an electron. The pion 
candidates in _B ^ ^07 must have Lk/{Lk + L^,) < 0.2 
and Lp/ {Lp + L„) < 0.5 and must not bo consistent with 
either an electron or a muon candidate hypothesis; in ad- 
dition, for all candidates with laboratory momenta above 
0.6 GeV/c, the number of photons observed in the DIRC 
is required to be consistent with the number that is ex- 
pected for the pion hypothesis. The performance of the 
pion identification requirements is evaluated with the de- 



cay _D*+ D'^{-^ K~tt^)tt~^ , which provides a large, 
clean sample of tt^ and . Using the results shown in 
Fig. 1 , we find the pion identification requirement retains 
85% of the pions from B ^ decays and rejects 99% 
of the kaons from B — > K*'^ decays. 

We form tt" candidates from pairs of photons with en- 
ergies greater than 50MeV in the laboratory frame and 
an invariant mass m^^ in the range 115 — 150 MeV/c^. 
We combine the identified pions into vector-meson can- 
didates requiring 630 < m^+„- < 960 MeV/c^ 640 < 
m^+^o < 930 MeV/c^, and 760 < m^+^-^o < 
790 MeV/c^ for p+, and o), respectively. The charged- 
pion pairs are required to originate from a common ver- 
tex. 

The photon and p/uj candidates arc combined to form 
the B meson candidates. We define /S.E = E^ — \/i/2, 
where E*^ = E*^^ + E* is the CM energy of the B meson 
candidate. The AE distributions of signal events are 
expected to peak near zero with a resolution of about 
50 MeV dominated by the photon energy resolution, and 
to have a tail in the negative region due to photon energy 
loss in the detector. We also define the beam-cnergy- 

substituted mass rriES = \J s/4: — p'^^, where p'^ is the 
CM momentum of the B candidate modified by scaling 
the photon momentum so that E*^^ + E'*^ — v^/2 = 
0. This procedure improves the tties resolution for the 
signal events in the AE negative tail. Signal events are 
expected to have an rriEs distribution centered at the 
mass of the B meson tob with a resolution of 3 MeV/c^. 
We consider candidates with tbes > 5.22 GeV /c? and 
—0.3 < AE < 0.3 GeV for further analysis. This region 
includes sidebands that allow the continuum background 
yields to be extracted from a fit to the data. 



B. BAGGED DECISION TREE 

The bagged decision trees are trained separately for 
the 5+ P^li -B" /3°7, and B^ uj-y channels with 
MC simulated signal and background event samples. The 
background sample consists of a BB MC sample that is 
about 3 times larger than the data and of a continuum 
MC sample is about 1.5 times larger. For the input clas- 
sifiers, we choose approximately sixty event quantities 
that characterize the kinematics of the candidates, 
the high-energy photon, the vector meson, the B me- 
son and the rest of event (ROE), which are the particles 
that are not used to reconstruct the B candidate. These 
quantities all have distributions that agree well between 
off-resonance data and continuum MC events. 

To reduce combinatorial background in the recon- 
structed TT^ candidates, we use in the BDT the invariant 
mass rrijj and cos6'-y-^, the cosine of the opening angle 
between the photons in the laboratory frame. 

We associate the high-energy photon candidate 7 with 
each of the other photons 7' in the event and calculate 
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FIG. 1: Performance of the charged-pion identification requirement applied to S — » p7 decays, evaluated using the D* control 
sample. Filled circles are for tt* efficiency and use the left-hand scale. Open circles arc for mis-identification and use the 
right-hand scale. The plot on the left shows results for continuum MC events and the plot on the right shows results for data. 



the likelihood ratio 

^ P(m-yy,.Ey|z) 

' Vim-yy, Ey \signal) +'P{m^Y,EY\i)' ^' 

where i = 7r°, r; and P is the probability density function 
(PDF) defined in terms of the energy of the second pho- 
ton in the laboratory frame Ej' and the invariant mass of 
the pair m^-y/. The PDFs are determined from simulated 
signal and continuum background events. The likelihood 
ratios CTZi^a and £7?.^ arc used in the BDT to reject high 
energy photons from 7r° and rj decays. 

To reject background events from B — > p(7r°/r?) and 
B —> uj{'!t'^ /r]), we also use the vector meson helicity an- 
gle, 9h, which is defined as the angle between the B 
momentum vector and the 77+ track calculated in the p 
rest frame for a p meson, or the angle between the B 
momentum vector and the normal to the lo decay plane 
for an uj meson. This variable is useful because in signal 
events, the vector meson is transversely polarized, while 
in the background events it is longitudinally polarized. 

Variables used in the BDT to reduce continuum back- 
ground include R2, the significance of the separation of 
the two B vertices along the beam axis (Saz), the polar 
angle of the B candidate momentum in the CM frame 
with respect to the beam axis (Og), and i?2, which is 
R2 in the frame recoiling against the photon momentum. 
We compute the moments Mj = J2jPj ' kos6'*|V Z^^ pj 
with i = 1, 2, 3, where p* is the momentum of each par- 
ticle j in the ROE and 6j is the angle of the momentum 
with respect to an ax;is. We use the Mj with respect to 
the photon direction and the ROE thrust axis. We also 
include flavor-tagging variables [11] to exploit the differ- 
ences in lepton and kaon production between background 
and B decays. 

While we find that all the variables contribute to the 
sensitivity of the analysis, the most effective ones are 



Saz, coaOy-,, R2, cosOg, M3 with respect to the pho- 
ton direction, the missing mass of the ROE, cosdn. and 
£TZ^o„. The distribution of the BDT output for the de- 
cay B" p°7 is shown in Fig. 2. We require the BDT 
output to be greater than 0.94 (0.93) for B ^ p-f (B° 
W7). These selection requirements have been optimized 
for maximum statistical signal significance with assumed 
signal branching fractions of 1.0x10^^ and 0.5x10^® for 
the charged and neutral modes, respectively. The signal 
significance is determined from a fit described in the next 
section. For the signal events that pass the loose selec- 
tion criteria, the BDT requirements have an efficiency of 
19% for B+ p+7, 31% for ^ p^, and 34% for 

In events where multiple candidates are present, we 
select the one with the reconstructed vector meson mass 
closest to the nominal mass. This criteria is chosen be- 
cause the mass of the vector meson is found to be uncorre- 
lated with the variables used in the fit. After applying all 
the selection criteria described above to signal MC sam- 
ples, we find signal efficiencies of 4.2% for B+ — * p"'"7, 
7.7% for B" p"7, and 5.2% for B" u}-f (taking 
into account the branching fraction B(u> 7r"'"7r~7r°) = 
0.892 ± 0.007 [12]), while backgrounds are reduced by 
0(10-5). 



IV. MAXIMUM LIKELIHOOD FIT 

We determine signal yields from an unbinned maxi- 
mum likelihood fit to toes and AE. The likelihood func- 
tion for a signal mode k (= p'^j, p°7, lv^) with a sample 
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FIG. 2: Distributions of the BDT output for B° p°j in 
signal (dashed) and background (dotted) MC samples. The 
distributions are normalized to 1. 



of Nk events is defined as 

/ ^hyp 

£fc = exp - ^ ; 



j=l \ z=l 



, (3) 



where A^hyp is the number of event hypotheses, and rii 
is the yield for each. For uj'j, three event hy- 

potheses are considered: signal, continuum background 
and combinatorial B backgrounds. For i?" p°7, a 
if*°7 background hypothesis is also included, 
while for _B+ P^l, a combined /f*+7/p+7r° 
hypothesis is included. Since the correlations between 
Toes and Ai? are found to be negligible in MC event sam- 
ples, we define the probability density function Vi{xj \ di) 
as the product of individual PDFs for each observable 
Xj=TOES, given the set of parameters on. 

The individual PDFs are determined from fits to ded- 
icated MC event samples. The signal toes PDFs are 
parametrized by a Crystal Ball (CB) function [13] and 
the /S.E PDFs are parametrized as 



f{x) = exp 



(4) 



where fi is the peak position of the distribution, ctl and 
aR are the widths on the left and right of the peak, and 
aL and an are a measure of the tails on the left and 
right of the peak, respectively. The peak positions and 
widths of the signal toes and AE PDFs are corrected 
for the observed difference between data and MC sam- 
ples of i? ^ K*^ decays. The PDFs for the remaining 
B° ^ K*°j and combined B+ K*+-//B+ p+n° 
backgrounds are determined from dedicated MC samples 
that are 100 times larger than the data. These PDFs 
are described by a CB function for toes, with a peak 



position the same as that of the signal PDF but having 
a much larger width, and a CB function for AS, with 
a peak position near —80 MeV. The negative AE peak 
position is due either to a kaon misidentified as a pion in 
B K*j or to a single missing photon in 5+ p'^ir'^. 
The Toes and AE PDFs for all other B backgrounds are 
determined from the BB MC sample. The toes spec- 
tra peak slightly in the signal region, and therefore are 
parametrized by a CB function, while the AE spectra 
are parametrized by an exponential function. The con- 
tinuum Toes and AE PDFs are parametrized by an AR- 
GUS threshold function [14] and a first order polynomial, 
respectively. 

The fit to the data determines the signal yield rigig, the 
continuum yield and the shape parameters of the con- 
tinuum PDFs. The shape parameters of the signal and 
B background PDFs are fixed in the fit. The relative 
yield between the peaking and the other B backgrounds 
is fixed to the value obtained from known branching frac- 
tions [12] and selection efficiencies determined from MC 
event samples. The overall yields of the B backgrounds 
are also fixed. All fixed parameters are later varied to 
evaluate systematic errors in risig. 

We validate the fitting procedure using ensembles of 
signal and background events. Two types of ensembles 
are produced: both signal and background events gener- 
ated using the PDFs described above; signal events ran- 
domly sampled from the GEANT4 MC events and back- 
ground events generated using the corresponding PDFs. 
No bias is found in the fit to these event samples. 

Figure 3 shows the data points and the projections of 
the fit results for AE and toes separately for each decay 
mode. The signal yields are reported in Table II. The 
significance is computed as v2AIn£, where A In £ is the 
log-likelihood difference between the best fit and the null- 
signal hypothesis. To take into account the systematic 
error in risig, the likelihood function is convolved with 
a Gaussian distribution that has a width equal to the 
systematic error. 



V. SYSTEMATIC UNCERTAINTIES 

Table III gives the contributions to the systematic un- 
certainties. The systematic error affecting on the sig- 
nal efficiency includes uncertainties on tracking, particle 
identification, 7 and 7r° reconstruction, and BDT selec- 
tion. The modeling of signal and background in the fit 
contributes to the uncertainties on the signal yields. 

The errors in BDT selections are determined from a 
control sample of the decay i?" K*^{-^ K^tt^)j for 
the mode and a sample of K*~^{-^ K^tt'^)j 

for the p+ and u. These B — > K*j decays are kine- 
matically similar to i? — > (p/uj)^ decays. The events 
are required to pass all applicable loose selection crite- 
ria, except the pion identification requirements. We also 
require the invariant mass 0.80 < to^+tt- < 1-0 GeV/c^ 
and 0.82 < TO^f+^o < 0.96 GeV/c^. The BDT output 
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FIG. 3: AiJ and rngs projections of the fits for the decay modes — > p"'"7 (top), p^j (middle), and — * uj^ 

(bottom). For illustrative purpose only, these plots are made by requiring —0.2 < AE < 0.1 GeV for the thes projections and 
"^ES > 5.27 GeV/c^ for the AE projections. The points are data, the solid line is the total PDF, the dashed line is the sum of 
B background PDFs, the dash-dotted line is the continuum PDF, and the dotted line is the signal PDF. 



classifiers are computed from the decision trees trained 
for the corresponding signal modes. The differences in 
the BDT selection efficiencies between the B — > K*'^ 
data and MC samples are used to correct the signal ef- 
ficiencies. The efficiency correction factor is 0.88 ± 0.09 
for B+ p+7, 0.91 ± 0.04 for B° p°-f and 0.90 ± 0.05 
for B'^ — !■ ujj. The uncertainty of the correction is taken 
as the systematic error. The large BDT systematic er- 
ror for the decay B^ /9+7 is due to the limited size 
of the 5+ — > K*^{-^ A'+7r'^)7 sample. As a means of 
validating the BDT technique, we apply the same anal- 
ysis technique to the B K*^ data control samples 
and measure the branching fractions for B K*"/. The 



results are consistent with the world averages [12]. 

The error in the pion identification requirements is es- 
timated using the D* control sample as shown in Fig. 1. 
Based on the difference of a momentum-weighted effi- 
ciency between the continuum MC sample and data, a 
1% systematic error per charged-pion is assigned to the 
B ^ decays. The MC sample is in better agree- 
ment with data for the looser pion identification criteria 
applied to i?*^ — > W7 and a 0.5% error per charge-pion 
is assigned. The uncertainties from tracking, tt" recon- 
struction, and photon selection are also determined from 
suitable independent data control samples. 

To estimate the uncertainty related to the modeling 
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TABLE II: The signal yield risig, significance E in standard 
deviations including the systematic error in rieig, efficiency 
e, and branching fraction B for each mode. The first error 

is statistical and the second is systematic. The branching 
fractions for B (p/iu)j and B pj are obtaintxl with the 
assuiuptioii of isospiii and Sr(-i)K syimiictrit-s. 



Mode 



€(%) 



B(10- 



B+ 
B° 



p+-i 23.3l7:3 ± 3.1 3.2(7 4.2 l-2Qtl fj ± 0.20 
34.9l^:^±1.2 5.4(7 7.7 Q.mtl%l±Q-Q& 



P°7 



B° uj-i 12.4l^;^ 


±2.0 


2.2(7 


5.2 


o.5oie;-^^ 


±0.09 


B ^ 




()..") (T 






± 0.16 






6.0(7 




1 70+0.34 

'^-0.32 


±0.17 


TABLE III: Fractional 


systematic errors (in %) of the mea- 


sured branching fractions. 










Source of error 


p+7 


p"7 


017 


P7 


(p/u;)7 


Tracking efficiency 


0.4 


0.4 


0.4 


0.4 


0.4 


Particle identification 


1.0 


2.0 


1.0 


1.4 


1.2 


Photon selection 


2.8 


2.8 


2.8 


2.8 


2.8 


7r° reconstruction 


3.0 




3.0 


1.7 


2.0 


BDT efficiency 


9.3 


4.2 


5.1 


7.0 


7.5 


Signal model 


7.1 


2.1 


16. 


3 3.0 


3.0 


Background model 


10.9 


2.8 


2.7 


4.3 


3.6 


BB counting 


1.1 


1.1 


1.1 


1.1 


1.1 


B{u> — > 7r"'"7r~7r°) 






0.8 




0.1 


Sum in <iua<lratur(> 


i(i.r 


().() 


17. 


) !).") 


9.5 



using a Bayesian technique, assuming a prior flat in the 
branching fraction and taking into account the system- 
atic uncertainty. 

We test the hypothesis of isospin symmetry by mea- 
suring the quantity 



r(E+ ^ p+7) 
2r(B0 ^ PO7) 



-o.43t^:2i±o.io. 



Most theoretical calculations [1-3, 15] predict small Ap. 
For example the estimate in Ref. [3] is —0.05 ± 0.03 for 

7 = 60° and -0.10±0.02 for 7 = 70°, where 7 is the phase 
of V*^. Our result is consistent with these predictions 
within the large experimental errors. However, it is worth 
noting that a recent calculation [16] indicates that non- 
perturbative charming penguin contributions can accom- 
modate large A^. We also measure the SU(3)F-violating 
quantity 



-0.49t^;|^±0.10, 



which is consistent with the theoretical calculations. 

We extract average branching fractions using a simul- 
taneous fit to all the relevant decay modes with the con- 
straints on the widths of the decay modes: V g^^pJr^^ = 
The average branching fractions 



2r, 



are defined as 



2r 



of the signal and background, we vary the parameters 
of the PDFs that are fixed in the fit within their er- 
rors. We vary the relative and absolute normalizations 
of B background components that arc fixed in the fit 
based on a kaon mis-identification study using the D* 
control sample as shown in Fig. 1. We find the difference 
in the momentum-weighted kaon mis-identification rates 
between the data and MC samples is 23% and conserva- 
tively vary the B — > K*^ background yield by 30%. The 
effect of the uncertainty of B{B^ P^""") [12] is also 
considered for the decay B^ p^7. For all the vari- 
ations, the corresponding changes in the extracted sig- 
nal yield are taken as systematic uncertainties, which are 
then combined, taking into account correlations. The er- 
ror on background modeling for 5+ p+7 is dominated 
by uncertainties in B background PDFs. 



B{B ^frf) = l \b{B+ ^ p+7) ± 2^B(B0 ^ 

(5) 

and 

B[B^{p/u;)^]^^[b{B+^P+^) 

where tb+ / tbo is the measured ratio between the charged 
and neutral B meson lifetimes, for which the current 
world average is 1.071 ± 0.009 [12]. Our measurements 
of the individual branching fractions are consistent with 
this hypothesis, with a of 2.3 for 2 degrees of freedom. 
We find: 



VI. RESULTS 

To calculate the branching fractions from the mea- 
sured signal yields, we assume B{T{4:S) B'^B^) = 
B{r{4S) B+B-) = 0.5. The results are listed in Ta- 
ble II. For B^ Luj, we also compute the 90% confi- 
dence level (C.L.) upper limit B{B° ^7) < 0.9 x 10-'^ 



B{B ^ p7) - {l.Tit'o.lt ± 0.17) X 10-6 
B[B^ip/cuh] = (1.63+°;30± 0.16) X 10-6. 

Using the world average value of — > K*~^^) = 
(4.03 ± 0.26) X 10-^ B{B" -> K*'^-f) = (4.01 ± 0.2) x 
10-^ [12], and the isospin averaged branching fraction 
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B{B K*j) = (4.16 ± 0.17) X 10-^ we calculate 



Rp+ 



RpO 



Rn 



■p/u 



B(B+ 


P+7) 

r 1 / 


B(B+ 


A'*+7) 




r I ) 


B(BO - 








B{BO - 




B{B- 


Pi) 


B{B^ 


■ K*j) 


B[B^ 




B{B- 





= 0.030 



-0.012 
-0.011 



= 0.024 ± 0.006 



U.U1Z_Q QQg 



0.042 ±0.009 



0.039 ± 0.008. 



isospin symmetry by 
= -0.43+0:2^^ ±' 0.10. 
branching frac- 
(l-73+o;32 ± 0-17) X 10-f' and 



We test the hypothesis of 
measuring the quantity Ap = 
We also measure the averaged 
tions B{B p-f) = (1.73t°-3'* 

B[B {pIoj)i] = (1.63t^J^ ± 0.16) x lO^^. These 
results are in good agreement with, and supersede, the 
previous published BaBAR measurement [6], which uses 
a subsample of the data used for this analysis. These 
results are also consistent with the measurements from 
Belle [18]. These branching fraction measurements are 
used to extract |Vtd/Vts| in a way that is complementary 
to the approach using B mixing [4]. 



These ratios of branching fractions can be used to calcu- 
late \Vtd/Vts\ [3, 5, 17]. Following Eq. 1 and using 1/Cp = 
1.17 ± 0.09, l/Co; = 1.30 ± 0.10 [3], Ai?p+ = 0.057^°;°^^, 
Aiipo = 0.006tHl3' and Ai?„ = -0.002t° °^^ [1] , we 
obtain 



\Vtd/Vts\p+ 

\Vtd/Vts\pO 

\Vtd/Vts\^ 



= 0.198^^35 ± 0-016 
= 0.254l^;°i^? ± 0.021 
= 0.202lg;^|g ± 0.016, 



where the first error is experimental and the second is 
theoretical. Using the average branching fractions and 
following Ref. [1], we obtain 



\Vtd/Vts\p 

\Vtd/Vts\p/ui 



= 0.235t!];^^^ ± 0.020 

— n oqq+0. 025+0. 022 



Similar values are found following Ref. [3]. These re- 
sults are consistent with the value of this ratio, 0.208 ± 

0.002(cxp)l°;°og(theory) [12], obtained from the studies 
of Bd and Bg mixing by the CDF and DO Collaborations. 



VII. SUMMARY 



We report the updated measurements of the branching 



fractions for the radiative decays 
and B^ — > uj^ 



P+7, S° 



B{B+ ^ p+7) = (1.20tgi^ ± 0-20) X lO'^^ 
^ p°7) = (0.97tg j^ ± 0.06) x lO'^^ 
B{B° w^) < 0.9 X 10-<5 (90% C.L.). 
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